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Abstract 

The Anti-de Sitter Space/Conformal Field Theory (AdS/CFT) correspondence may offer 
new and useful insights into the nonperturbative regime of strongly coupled gauge theories 
such as quantum chromodynamics (QCD). Recently a modified soft-wall AdS/QCD model 
incorporated independent sources for explicit and spontaneous chiral symmetry breaking and 
linear confinement. This model contains a modified dilaton and higher-order interaction terms 
in the Lagrangian. Within this model we explore the radial pseudoscalar mass spectrum using 
two different representations of the pion field. We find the mass eigenvalues associated with 
each representation, show the equivalence between the two, and find good agreement with the 
pion masses. The Gell-Mann-Oakes-Renner (GOR) relation is naturally obtained. 
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1 Introduction 



The theoretical framework for Anti-de Sitter Space/Conformal Field Theory (AdS/CFT) corre- 
spondence, relating type IIB string theory in AdS$ x S 5 to J\f = 4 Super Yang-Mills (SYM) theory, 
was laid out in (TJ-Sl. This correspondence led to an effective dictionary relating strongly coupled 
gauge field theories to higher-dimensional, weakly coupled gravity theories. Such a correspon- 
dence removes the main difficulty of strongly coupled systems, the inapplicability of perturbation 
theory. Many papers have constructed dual models to capture part of the nonperturbative regime 
of gauge theories [fSMTOll . The majority of these use the correspondence dictionary to describe 
quantum chromodynamics (QCD)-like theories containing a large number of colors N c . 

Two methods are employed to formulate models using AdS/CFT correspondence: top-down 
and bottom-up. A top-down method uses a version of string theory to calculate an effective La- 
grangian that, it is hoped, will contain certain key characteristics of QCD. A bottom-up approach, 
commonly known as AdS/QCD, uses the basic tenets of QCD in d dimensions to formulate a 
dual gravity theory in AdS d+ i. Of course, constructing an AdS dual theory that encompasses the 
richness of QCD presents the greatest challenge, a task yet to be accomplished. In this paper, 
we implement a bottom-up approach evolving from the work of ll5l-l7l[TTTl. often associated with 
the term holographic QCD. This phenomenological model incorporates many of the crucial as- 
pects of QCD. While the exact gravity dual is not known at present, nor even whether one exists, 
these types of models spark interest because of the potential of uncovering universal properties of 
strongly coupled gauge theories. 

We continue to study the model presented in IfTTTl by exploring the pseudoscalar sector. That 
model improves upon the previous ones both by incorporating confinement and separating explicit 
from spontaneous chiral symmetry breaking as in QCD. We are particularly intrigued by this prob- 
lem because the pseudoscalar representation used in (6l[l2l appears to be incompatible with the 
representations used in ElMIMIlll. Both are in the exponential form 

X = Se iP , (1) 

where P is the pseudoscalar part of a complex field X. The former papers define P using the 
vacuum expectation value v(z) of the scalar field, S(x, z), as P = n a (x, z)t a /v(z), while the latter 
papers specify P = 2n a (x, z)t a . (The x = (t, x) while z is the fifth-dimension coordinate.) The 
two representations produce seemingly different equations of motion and potentially different mass 
eigenvalues. Alternatively, one may use a linear representation, X = X + i7r a t a , as in lfi~5ll . Any 
physical observables calculated from this model should of course be independent of the choice of 
representation. 

In this paper, we begin by recalling the modified gauge-gravity dual model in Section [21 For 
illustrative purposes, we examine the most prevalent exponential representation and then compare 
it to a linear representation. In Section [3] we derive and explore the equations of motion of the 
two representations, showing that the sets of equations are, indeed, equivalent. One can simplify 
the equations of motion in two different ways: (i) by eliminating one field as in |fi~6l , or (ii) by 
transforming the full system of second-order differential equations into Schrodinger-like form. 
Numerically calculating the mass eigenvalues in Section we find that the eigenvalues found 
from one method do not match those of the other. Apparently, lfi~6l has found the eigenvalues 
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not of the pseudoscalars but of a derivative field, and they are not equivalent. The pseudoscalar 
mass spectrum calculated from the model of fl~Q agrees reasonably well with the observations. In 
Section [5] we verify quantitatively that the pion mass-squared varies linearly with the quark mass, 
thus satisfying the Gell-Man-Oakes-Renner (GOR) relation. We conclude with general comments 
on the pseudoscalar representation in Section [6l 



2 The Dual Model 



We investigate a modified version of a common soft-wall AdS/QCD model. The original model 
was introduced in and further investigated in ll8 UT0l[T3l[T7l[T8l . A modified dilaton profile and 
quartic interactions were introduced in ifTTI to separate explicit from spontaneous chiral symmetry 
breaking. The action for this AdS/QCD model is 



S 5 = - j d 5 x^e~ x{z) Tr 
where the metric takes the form of 



\DX\ 2 + m 2 x \X\ 2 



k\X\ 4 + 1 ±- 2 (F 2 a + F } 

z !J5 



ds 2 = a[z) 2 {ri^dx^ dx v + dz 2 ) 



(2) 



(3) 



with the ( — h + + +) signature and the warp factor a(z) = L/z with z > 0. The covariant 
derivative is defined as Dm = 8m + i\Vj\t, X] — i{Am, X} (upper-case roman indices run from 
to 4). The complex 2x2 field X contains the scalar vacuum expectation value (VEV) v(z), 
scalar field excitations S(x, z), and the pseudoscalar field n(x, z). Using the AdS/CFT dictionary, 
the scalar component is dual to the qq operator and, therefore, the scalar mass associated with this 
field in the 5-D gravity theory is m x = — 3/L 2 . We are only interested in a subset of terms related 
to the pseudoscalar sector, 



\DX\ 



m 2 x \X\ 2 + K\X\ 



2g\ 



:{d M A N - d N A 



M i 



ge -x(*) Tr g MN (d M X - i{A M ,X} - i[V M ,X])(d N X* + i{A N ,X*} + i[V N ,X^)) 



+m x \X\ 2 - k|X| 4 + 



gMP 



(d M A N d P A R - d M A N d R A P ) 



A good phenomenological parameterization was described and justified in IfTTI . 



v{z) 



az + /?2;tanh(7z 2 ) 

d z {a{zfv\z)) + a{zf (^kv(z) 3 + 3v{z)) 



dz ■ 



a(z) 3 v'(z) 



(4) 

(5) 
(6) 
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where 



a 




(7) 



95L 




V kL 2 

95° 



4A 



— a 



(8) 



7 



(9) 



For the rest of this paper we use the values quoted in [fTT|| : m q = 9.75 MeV, a = (204.5 MeV) 3 , 
A = 0.18 GeV 2 , g 5 = 2tt, and k = 15. 

Concentrating on the pseudoscalar sector, we explore two representations for the field X, 



where / is the 2x2 identity matrix. We refer to X e as the exponential representation and Xi as the 
linear representation. The exponential representation is used in (51[8l[T3l[T4l, where it is assumed 
to be the canonically normalized pion field, n e = it/ f^. The representation X e is also used in lfT6l . 
where the eigenvalues of the pseudoscalar sector are computed using a method that we comment 
on later. The linear representation has been specified before in [fT5l . where the tt field carries the 
same dimensions as other fields in the Lagrangian. We already see an apparent difference between 
representations; (fTTI) allows for an explicit quartic term in n when substituted into the action ©, 
whereas there is no such term in the case of (flOi The consequence of such quartic terms in n will 
not be addressed here; however, the quartic term strength impacts the pseudoscalar mass spectrum 
through the parameter k. We only consider field terms up to quadratic order. 

3 Representations 

In this section we derive the equations of motion arising from the two representations (flOl) and (fTTI) . 
The pseudoscalar and longitudinal components of the axial-vector field mix in the Lagrangian; 
therefore, we find two coupled differential equations for each representation. This makes the 
numerical work more involved than for the scalar, vector, and axial-vector sections, which were 
already studied in [fTTI . The last part of this section shows that the two sets of differential equations 
are equivalent. 

3.1 Exponential Representation 

Let us take (fTOl) and substitute it into ©, where we focus only on the terms involving the field 

ir(x,z), 



X e 
X,. 



(v(z)/2 + S(x, z)) Je 2<f °t I# 
(v(z)/2 + S(x, z))I + irf(x, z)f 



a 



(10) 

(11) 





(12) 
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We work in the axial gauge, A z = 0, and define = A pl _ + <9 M 0, where d^A^ = 0. Separating 
(PT21) explicitly into regular 4-D components and extra-dimensional terms, we obtain 



-qq^(v 2 d u nd u 7i + v 2 A u A u - 2v 2 d u 7tA l 



gg zz v 2 d z 7rd z 7r 



9l 



(d^A p d u A a - d^ApdvAv) + 



-gg zz g^ 
gi 



{d z A»d z A u ) 



Keeping only terms of the longitudinal part of A p gives 



— e~ x{z) 



-gg h 



^(v 2 d^d u Tr + v 2 d, 



2v 2 d^7id l/ (, 



+V^gg zz v 2 d z nd z 7r + 
Varying (fT4l) with respect to it gives 



-gg zz g» u 
gl 



5C e = d z e- x ^=g~g zz v 2 d z -K5>K + e^V^V^^TT - 0)<5vr. 
Using a Kaluza-Klein decomposition, 

ir(x,z) = y^n w (x)7r ra (^) 

n 

<fr(x,z) = ^2<& n (x)(j) n (z) 



(13) 



(14) 



(15) 
(16) 



and 



d 2 U n (x) = m 2 n U n (x) , d 2 $ n {x) = m 2 J> n {x) 



(17) 



we can express the system of equations in terms of its 2-dependent parts. We obtain the first 
equation of motion, 

^—d z (e- x v 2 a 3 d z rr n ) + m 2 n (rr n - <j> n ) = 0. (18) 

Varying (fT4l) with respect to and breaking it down into KK modes gives the second equation of 
motion, 



e x d z 



z 



-d z <p n + 



7T„ - 



n fn J 



0. 



(19) 



Alternatively, we can express (1181) and (1191 ) in a Schrodinger-like form. We rewrite them in the 
same form as (|64|) and (1651) in the Appendix by substituting 



Tin 



e f{z) n r 



g{z) = x( z ) + ^gz, 



(20) 
(21) 
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which eliminate terms involving the first derivative of the fields, ir' n and (j)' n . After reverting back 
to the notation tt — > ir and — > <p the equations of motion become 

„ fx' 2 x" X'v' 3x' 15 3v' v" 2 \ 2 v 2 L 2 ± ^ 

_7F n +\— A o + 1T + 71 + m n)^n = -7Tl n ^—(j) n (22) 



4 2 v 2z 4z 2 vz v n J n z 2 

" l x ' 2 x \ x ' 3 + ^ W> U o 2 rr nx\ 
+ ( T - T + ^ + ^ + ^ 2 -J^ = ^ (23) 

where (') indicates the derivative with respect to z. 



3.2 Linear Representation 

When considering the linear representation of the pseudoscalar field (fTTI) . we find quadratic and 
quartic n terms that were not explicitly present in the exponential representation. After making the 
appropriate substitutions, we find that 

+ g^v 2 d^d v <P+ y —^-d z d^d z d u <py (24) 

Once again, we derive two coupled equations. Varying with respect to produces a result similar 
to X e with the exception of factors of the VEV in the mixing term, giving 

e x dz l—d z <f> n ) + (vr n - v<f> n ) = 0. (25) 

Varying with respect to n gives the second equation of the linear representation, 

Z 3 e X O z ^^-<9 Z 7T„^ - - —T^T^j n n + m n 7r " = m n^^' ( 26 ) 

We can express (|25T) and (|26|) in a Schrodinger-like form with the substitutions, 

T\ n = e / vf n / = - + - log - (27) 

K = e>; ^ = | + ilog£ (28) 

Simplifying the equations and reverting back to the notation ix n — > ir n and 4> n — > n for simplicity, 
we find 

-^ + (T-T + i? + l + ^-J^ = V r - (29) 

+ A + = 00) 



4 2 ' Az 2 ' 2z 2z 2 n n ~ n z Yn 



5 



3.3 Representation Equivalence 

The pseudoscalar field representation should not affect the physical results obtained from the 
model. Examining the two sets of coupled equations in each representation, we see that neither 
(TT8l) nor (fT9l) contains an explicit dependence on k, whereas (1261) does. Although dependence on 
K does not appear explicitly in the exponential representation, it does enter through the function 

v{z). 

We begin by expanding X e , 

X e = (| + 5)(l + 2i7r e + ...) 

= V - + S + m e v. (31) 

Comparing (|3TT) to (fTT|) . we surmise that n e v(z) — > t\\ relates the two representations. Let us 
substitute ix e — > iti /v(z) into the equations of motion of the exponential representation and attempt 
to obtain the equations of motion of the linear representation. The substitution into (fl9l ) is trivial; 
it yields 

e x dz (^d z( yj + ^fa - v<j>) = 0, (32) 

which is equivalent to (|25T ) as expected. Showing the equivalence of the other two equations 
requires a bit more work. First, we substitute for ix e in (fT8l) and then simplify the expression, 

z 3 e* fe~ x v 2 f-n' TXiv\\ 2 , 

— a ' {— mUni ~ V(p) = (33) 

which becomes 

< ~ (x' + < - ^ (v" - X'v' - + ml^ - v<j>) = 0. (34) 

Recall the equation of motion for v (z), found in [1 1], which can be derived from (0 and does not 
depend on the representation. 

„ / , 3\ , /3 kL 2 v 2 \ 



Using d35j) in C34J, we find 



<-\z + x )<+[- + )*i + K fa - v(f>) = 0, (36) 



3 A , /3 kLV 

which is the same as the equation of motion of the linear representation (|26l) . In a similar way, this 
equivalence can be shown by starting with the linear representation and substituting ni = v(z)-K e . 
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4 Pseudoscalar Mass Eigenvalues 



We investigate two ways to calculate the pseudoscalar eigenvalues rr? n . Rearranging and elimi- 
nating the longitudinal component <p is one strategy outlined in [fT6ll and is briefly presented here. 
Alternatively, the Appendix contains a numerical routine we use to solve the set of coupled equa- 
tions. Using this method, we find that 7r e , the ratio of i\\ and v(z), is extremely sensitive to boundary 
conditions, the reason being that v(z) goes to zero as z goes to zero. This makes resolution of the 
eigenvalues difficult and subject to significant numerical error. Fortunately, we have shown ex- 
plicitly that physical results do not depend on the particular representation of the pseudoscalars. 
Therefore we determine the eigenvalues in the linear representation. 

First, we follow the method of [fT6ll . We manipulate (fT"8l) and (IT9T ) to eliminate the </> field. 
Adding (Q3]> and (O yields 

g 2 a 2 v 2 d z n n = m 2 n d z (j) n . (37) 
We then use ([37]) to replace the first term in (TT9b and solve for 0, 



1 -8, 



- \ I 95 2 2 a 

ae x — -a v d z i\ n 
T 



+ 7T n 



glK z ) 

m 2 n d z (j) n = d z (h(zy 1 d z (h(z)d z ir n )) + m 2 n d z ir n (38) 



where h(z) = a(z) 3 v 2 e x . Using (1371) again, we denote <9 2 7r — > n and rearrange to find the 
eigenvalue equation 

- d^h^d^hnn)] + g 2 a 2 v 2 TT n = m 2 n Tr n . (39) 



This can be put into the Schrodinger-like form by substituting U = n / ^Jh(z). Then (l39l) becomes 

- + V(z)U n = m 2 n U n (40) 

where the potential takes the form, 

3h' 2 h" g 2 L 2 v 2 
V & = -W~2h + ^- 

3 3v' v' 2 3 X ' v'V X' 2 v" X " 9 2 L 2 v 2 
= + 2— + -^- — + - A + 9 ■ (41) 

Az z zv v z 2z v 4 v 2 z z 

Solving (|40l) using a standard shooting method gives the mass spectrum shown in Table [Q There 
is no low-mass Goldstone boson and no large mass gap between the first two eigenvalues. What 
we have done is taken 2 second-order differential equations and reduced them to 1 third-order 
differential equation and found the eigenvalues of d z w. We seem to lose the Goldstone boson using 
this method. 

The numerical routine described in the Appendix calculates the mass eigenvalues of (|29l and 
(|3Q|) . which are then plotted in Figure [Hand listed in Table [T] Solving the set of equations directly 
produces a mass spectrum with a Goldstone boson and with mass eigenstates that match well with 
the observed radial pion excitations. These results show that eliminating one of the fields from the 
system of second-order differential equations also eliminates information from the mass spectrum. 
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2 4 6 

n 



Figure 1 : The mass spectrum calculated in the AdS/QCD model is plotted along with the experi- 
mental data |fi9ll . The eigenvalues display two characteristics matching the QCD pion spectrum: 
(1) low-mass ground state and (2) a large gap between the ground state and the first excited state. 
The large-n mass trajectory clearly follows our calculated eigenvalues from n « 4 when our nu- 
merical routine inadequately follows the oscillations of the higher eigenfunctions. 



n 


7r Data (MeV) 


txi (MeV) 


Large-n tti 


d z ir e (MeV) 


1 


140 


143 




1440 


2 


1300 ± 100 


1557 




1706 


3 


1816 ±14 


1887 




1925 


4 


2070* 


2095 




2117 


5 


2360* 


2298 


2245 


2290 


6 






2403 


2451 


7 






2551 


2601 



Table 1: The observed masses [fT9l and calculated masses using the linear representations. The 
large-n limit solutions are valid from n ~ 4. From that point onward, the numerical method used 
becomes increasingly inaccurate and tends to skew the m eigenvalues to larger values than are 
expected from the linear behavior. The eigenvalues found using the method of [fT6ll are also shown. 
* Appears strictly in the further states of |fT9l . 
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Thus, the eigenvalues of d z ii e do not match those of tti or ir e , at least not in this model. The validity 
of equating the two sets of eigenvalues in |[T6ll may need to be reassessed. 

For large-n excitations, the technique described in the Appendix runs into boundary condition 
problems for n > 4. As the number of oscillations in the eigenfunctions increases for higher 
n modes, the routine finds eigenvalues that are skewed to larger values. To uncover the correct 
asymptotic behavior for large n, we take the large-2; limit of (1291 and (|3Qb . As n increases, the 
eigenf unction is largely determined by the behavior of the effective potential at large z. At large z, 
the dilaton and tachyon behave as 

v(z) = (a + /3)z = Ty (42) 
X(z) = A^ 2 . (43) 

To take the large-z limit of both representations, we introduce a new dimensionless parameter, 
£ = \/Xz, and expand in £. In the exponential representation, we find that ((22)) and (1231) become 

2 

777 

-< + ^k = -f(7i k -T4> k ) (44) 
-0fc + £ 2 0fc = 9 ^{ix k -Y<p k ). (45) 



Similarly, in the linear representation the expansion of (|29l) and (1301) at large £ yields the eigenvalue 
equations 

» 9 / ^r 2 rnl\ mlT 

+ = ("2A - 2 + x) 7rfc "~X~^ (46) 

-tt + e^ = ^K-r0 fc ), (47) 

where (') indicates differentiation with respect to £. Each set of equations appears to describe a 
pair of simple harmonic oscillators, the equations of motion of which are 

-0' fc ' + £ 2 0fc = (2k + 1)^ (48) 
-< + ^k = (2k+l)ir k k = 0,1, (49) 



It is a reasonable assumption that <p k = c k n k ; it ensures that (1441) . (1451) . (1461) . and (1471) have solu- 
tions. Using the form of (|48l) and (|49l) to solve for c k and m 2 in both representations, we find 

ml = (2fc + 1)A + 9 |r 2 (51) 
for the exponential representation and 

Ck = gir* +%k + i)A (52) 

((2fc + 3)A-^r 2 ) (^ 2 r 2 + (2fc+l)A) 
mfc = (2jfe + l)A (53) 



9 



for the linear representation. 

So far, we have neglected an important fact: z > 0. The eigenfunctions <p and n describe half- 
harmonic oscillators and contain only half the modes that full harmonic oscillators do; therefore, 
we must take k — > 2k. From ifTTfl we have 

r 2 = — . (54) 

K 

The mass eigenvalues for large n, where n — k + 1, in both representations then become 

m 2 n = (An - 3)A + g 2 V 2 n = 4, 5, . . . (55) 

which are also listed in Table \T\ and plotted in Figure \T\ Combining (1551 ) and the numerical tech- 
nique, we obtain all the pseudoscalar eigenvalues. By simple investigation, we find that the large-n 
eigenvalues should be trusted over the ones found with the numerical routine for n > 4. 



5 Gell-Mann-Oakes-Renner Relation 

In this section we explore the Gell-Mann-Oakes-Renner relation numerically. Using the estab- 
lished equivalence between the exponential and linear representations, n e = iri/v (z), and inserting 
it into (1371) . we obtain 

9lLV aJ^]= m ldA. (56) 



2 z \ I ""ir^z 
Z z \ V 1 



Following the method of [5], we construct a perturbative solution in where (f>(z) = A(0, z) — 1 
and use 



2 r d 2 A(0,z) 

Jit — ~ Ll 2 

From this, it follows that 



(57) 



2 / n f , u d z A(0,u) 
ir(z)=mlv(z) du z \ ' ' . (58) 

Jo v ( u ) 9EU 

The function u 3 /v 2 (u) is significant only at small values of u ~ ^Jm q /a, where we may use (1571) 
to relate the derivative on A(0, u) to the pion decay constant, so that 

B = _!^fi. (59) 
v 2m q cr 

We find that i\i = —v(z) solves the axial- vector field's equation of motion found in [fTTTl 

( e' x \ n 2 o 2 L 2 v 2 

e x d z —d z A,(q, z) - 'LA^q, z) - — o — A^q, z) = (60) 
V z I z z 6 



in the region of small z and as q — > 0. As a result, (1591) becomes the expected Gell-Mann-Oakes- 
Renner (GOR) relation, 

2m q a = m 2 J 2 . (61) 
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4 6 
m (MeV) 
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Figure 2: Plot of m\ versus m q produces a straight line from which the pion decay constant f n is 
calculated using (l63l) . 



This perturbative behavior for n e and 7T; justifies the use of Neumann and Dirichlet boundary 
conditions, respectively, such that 



7T e (0) 



7T Z (0) = V(0) = 



The ratio of m q and should be a constant following the GOR relation. 

J IT 



mi 



2a 



(62) 



(63) 



We solve the pair of coupled differential equations for the ground-state pseudoscalar mass, m^, 
for differing values of m q to ensure that the numerical routine of the Appendix respects the GOR 
relation. The results are plotted in Figure [2l We see linear behavior in the plot, indicating that as 
m q — > we obtain a constant ratio of m q jm\. The slope of the line in Figure [2] implies f n m 90 
MeV, a result consistent with the input parameters as described in IfTTTl . 



6 Discussion 



We investigated two possible representations of the pseudoscalar field within a model of AdS/QCD 
that incorporates both explicit and spontaneous chiral symmetry breaking as well as confinement. 
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We showed formally the equivalence of these representations. The Gell-Mann-Oakes-Renner 
relation was shown to hold both analytically and numerically. We also found that the eigenvalues 
of the derivative field d z ix e are not the same as those of the field ix e , as IfTBI suggests. 

This paper completes the calculation of the radial mass spectra of the spin-0 and spin-1 fields in 
the soft-wall model of ifTTTl . The pseudoscalar fits into the AdS/QCD model quite well, reproducing 
a number of pion features, namely, (1) a low-energy ground state corresponding to the Goldstone 
boson, (2) a relatively large gap between the ground state and the first excited state, (3) a linear 
trajectory, and (4) a constant ratio of m 9 /m^. that gives a value of f n within a few percent of its 
accepted value. 

Acknowledgments TK thanks Brian Batell and Todd Springer for insightful discussions. TK 
also thanks Daniel Sword for help with the numerical techniques used in this paper. This work was 
supported by the US Department of Energy (DOE) under Grant No. DE-FG02-87ER40328. 



A Numerical Routine 

The equations of motion can be reduced to a set of second-order differential equations 

-<j ) " + V 1 (z)<J ) + f(z)7r = 
-tt" + V 2 (z)n + g{z)<j> = 



(64) 
(65) 



where the eigenvalues are contained within the coefficient functions. These equations can be reex- 
pressed as a system of first-order differential equations 



$' + W(z)$ = 



where W is the matrix 



W 



and $ is the vector 



/ o 


1 





o\ 


V x {z) 





/(*) o 











1 


\ g(z) 





V 2 (z) J 




/ 










-ft 














v - 


-<) 





(66) 



(67) 



(68) 



that forms an orthonormal basis of solutions. We can propagate the solution $ between two bound- 
ary points 

$(z 1 ) = U(z,z 1 ,z ,ml)$(z ) (69) 

where we solve (l66l) with the appropriate boundary condition at z . The eigenvectors and eigen- 
values of U are then calculated. We find two large and two small eigenvalues corresponding to two 



12 



nonrenormalizable and two normalizable eigenfunctions, respectively. Let us assume the eigen- 
vectors % and w 4 correspond to the small eigenvalues, A 3 and A 4 . Then, any solution for $j can be 
written as 

$i = au 3 + (3u 4 , (70) 



where we take the boundary condition as 

$iOo) 



(71) 



U3 u\ \ I a 



In order for a and (3 to be nontrivial, we must satisfy 
for Dirichlet or 



id u\ J V P 



(72) 



ut < ;ur° (73) 

for Neumann boundary conditions. We do this by cycling through eigenvalues w? n that minimize 



the determinant of the 2x2 matrix in (1721) or (1731) . Practically, we find the singular points in 
the graph of the quantity u\u\ — u\u\ (or — M3M4) versus m 2 n . An abrupt change in its 
behavior signals an eigenvalue. Of course, the elements chosen from the eigenvectors U3 and M4 
are dependent upon the choice of Neumann or Dirichlet conditions on the boundary zq. 
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